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Abstract 

We show that the recently introduced L 1 TV functional can be used to ex- 
plicitly compute the flat norm for co-dimension one boundaries. Furthermore, 
using L 1 TV, we also obtain the flat norm decomposition. Conversely, using the 
flat norm as the precise generalization of L X TV functional, we obtain a method 
for denoising non-boundary or higher co-dimension sets. The flat norm decom- 
position of differences can made to depend on scale using the flat norm with 
scale which we define in direct analogy to the X^TV functional. We illustrate 
the results and implications with examples and figures. 

1 Introduction 

In this research announcement, we point out that the L 1 TV functional, introduced 
and studied in the continuous setting in [3, 14, 1, 15] and earlier in the discrete setting 
in [2, 10], provides a convenient way of computing both the value of, and the optimal 
decomposition required by, the flat norm from geometric measure theory. 

The L X TV functional was introduced as an improvement to the now classic Rudin- 
Osher-Fatemi total variation based denoising [1 1]. Among its many nice properties are 
the way it handles binary images, making it useful for shape processing. Theoretically 
speaking, the clean geometric structure of the functional and its minimizers is very 
attractive. See [14, 1] for details. 
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Joan Glaunes was, as far as we know, the first to suggest and use the flat norm 
from geometric measure theory as a distance in shape space. In his dissertation [ 
and a couple of application papers [13, 6] with collaborators, the dual formulation of 
the flat norm is used to compute distances between shapes. 

Acknowledgements: It is a pleasure to acknowledge Bill Allard and Bob Hardt 
for useful discussions and Joan Glaunes and Sarang Joshi for inspiring this work. 
Additionally, the second author acknowledges useful discussions with Selim Esedoglu. 

2 L lr TV gives the Flat Norm 

The L 1 TV functional introduced and studied in [3] is given by: 

F(u) = / \Vu\dx + \ / \u — u \dx (1) 

where u and u are functions from IR n to K with u being the input function. In 
image analysis applications, n is usually 2 and u is the measured image intensity 
function. The optimal u minimizing (1) can be thought of as a denoised version of 
Uq. Typically, one chooses the parameter A based on noise levels since this choice 
selects the scale below which features or oscillations are ignored. In the case that the 
input function is binary, Chan and Esedoglu observe that the functional reduces to: 

F^ E (£) =Per(£) + A|£ Afi|. (2) 

where the binary minimizer is Xs ; the binary input is xn, Per(S) is the perimeter of 
S, and A denotes the symmetric difference. Of course xe denotes the characteristic 
function on E, with a value of 1 on I? and on the complement of E. Now, let 
A) be a binary minimizer of (2). I.e. 

E(fi, A) = argminF^(S) = Per(S) + A|S A Vt\. (3) 

For our convenience, we record the optimal decomposition of Q into A) and 

(£(fi, A) A fi)} as the pair {dQ, £(fi, A) A fi}. 

In what follows, we use the notions of current, mass and d and supporting ideas 
from geometric measure theory. We introduce and explain these in some detail in 
the Appendix. Informally, one can gain much by thinking of the n- current Te as 



2 



an orientable n-submanifold or n-rectifiable set E C M. n+k with an orientation, of 
the mass M(T#) as the n-dimensional volume of E, and of dT, the boundary of T, 
as the oriented boundary of E with orientation imposed by the orientation of E. 
We sometimes omit the subscript indicating the support of the current, referring to 
currents T and S. For those without experience with currents, we suggest focusing 
the examples in the the appendix. 

The flat norm of an n-current T, denoted by F(T), is given by 

F(T) = min{M( 1 S) + M(T — OS)}) (4) 

S 

where S varies over n + 1-currents and M is the mass of the indicated currents. We 
refer to {T, S} as the flat norm induced, optimal decomposition. Now for results. 

Theorem 2.1. For the current Tqq, 

¥(T dn ) = Fh B p(n, 1)) (5) 

and 

{Tgn, »?s(n,i)Aa} (6) 
is the optimal decomposition required by the flat norm. 

Theorem 2.1 says that L1TV computes the flat norm. This relation between the flat 
norm and the L1TV functional immediately suggests a very useful generalization of 
the flat norm. 

Definition 2.2 (Flat Norm With Scale). 

F\(T) = min{AM(5') + M(T — dS)}) (7) 

Theorem 2.1 is then simply a special case of the next theorem. 
Theorem 2.3. For the current Tga, 

W x (T dQ ) = F£ E (E(Q,\)) (8) 

and 

{Tan, Ss(n,A)An} (9) 
is the optimal decomposition that the flat norm with scale requires. 
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Proof of Theorem 2.3. The proof is really simply checking that the picture one 
can draw holds after the definitions of mass (M) and perimeter (Per) are used to 
translate the picture into analytic terms. Very briefly, we have: 



A|S A n\ +Per(S) 
XM(S^ n ) + M(T 9S ) 



Figure 1 illustrates this pictorially. 




□ 




Tan = Tqq - <9SV2ae 

Fc E (Y>) = A|E A n\ +Per(E) 



= \M(S SAn ) + M(T SS ) 

= \M(SsAto) + M(Tan - dS SA n) 



Figure 1: In this figure we illustrate the translation of the L X TV view to the flat 
norm view. The perimeter of S becomes the mass of Tqq — dS-£ A n and the volume 
of E A Q becomes the mass of S^aq- Note: this figure does not depict a minimizer. 
Rather, we depict Q and any candidate S. 
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Our final observation is that the flat norm with scale gives us the same decom- 
position as we would get if we first scaled T, computed the flat norm decomposition 
and then reversed the scaling. More precisely 

Lemma 2.4. Denote the optimal¥\ decomposition by {T, S}\. Then 

{T,S} x = d {# {d x# (T),d x# (S)}i (10) 
where d\ denotes the X-dilation ofM. m , and d\#(TM) = T^m)- 

Proof. The minimizing decomposition {T, S}\ which minimizes AM(S') +M(T — dS) 
also minimizes \ n ~M(S) + A n_1 M(T — dS) = M(d x #S) +M(d x# (T - OS)). Therefore, 
to get the minimizer for Fa run the optimization required by Fi using d\#(T) as input, 
and then contract with din. □ 

We now discuss applications and examples with pictures which should clarify 
things for those with less exposure to geometric measure theory. 

3 Applications and Illustrating Examples 

The value of the above results is fully realized by exploring their use in applications 
to images and shapes. As observed by Glaunes and others, the flat norm is a very 
natural candidate for distances between shapes. We now very briefly explore and 
illustrate applications of the above observations. 

3.1 Generalized Flat Norm: flat norm with scale 

As mentioned in section 2, by letting the A ^ 1, one has a natural way to vary the 
intrinsic scale in the flat norm. 

¥ X (T 9E ) = mm{XM(S)+M(T dE -dS)} 
= F£ E (E) 

where Tqe represents dE. The main point here is that this is easy to compute, given 
the connection to the L X TV functional. Varying A gives us the ability to choose what 
scale is big and worth keeping. See Figure 2. 
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JIM JUJOL 

A > 1 A < 1 

Figure 2: Flat norm with scale decompositions as a function of A: The larger A is, the 
less smoothing there is. Let T denote the blue input 1-currents and S denote the red 
optimal 2-currents. Then in fact, A is the bound on the curvature allowed in T — dS, 
see [1, 14] for details. Accordingly, the curves T — dS are allowed greater curvature 
in the curves on the left than is permitted in those on the right. 

3.2 Flat norm via L lr TV 

We can use the L X TV functional to very easily calculate both the flat norm of differ- 
ences between surfaces which are boundaries and the optimal decomposition of that 
difference into surface and area parts. In Figure 3, the decomposition of a boundary 
current Tqq into the diminished boundary Tqy, and the symmetric difference current 
Sj:m is illustrated. 

3.3 L lr TV by the dual form of the flat norm 

The dual formulation of the flat norm can be used to compute L X TV minimizers. In 
what follows, we define spt T to be the support of T. The following results establish 
that maximizing forms, or in some cases, maximizing sequences of forms, contain the 
decomposition into S and T — dS. 

Proposition 3.1. Suppose that T(<f>) = F(T) = min{M(S) + M(T - dS)), where 
(J) is a smooth, compactly supported n-form satisfying \(f>\ < 1, \d<f>\ < 1, and T is an 
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T = T dn 
and 



T; 



T-dS 



Figure 3: A simple reinterpretation of the L l TV input and minimizer gives us the 
flat norm of Tgn and the decomposition into the diminished boundary Tq% and the 
symmetric difference current S-£,txi- 



n-rectifiable current with density = 1 7i n almost everywhere on spt T, then on the 
support of S, \d(f)\ = 1, and on the support of T — dS, \<p\ = 1. 

Proof. For a minimizing choice of S, we have that 

T(<f>) = F(T) = M(S) + M(T - dS) 
T(4) = 6S(4>) + (T - dS)(<l>) 

and 

ds(cf)) = s(d<f)) 

so that 

S(d<j>) + (T - dS) (0) = M(S) + M(T - OS) (11) 

We know M(S) = J ldH n+1 LsptS and M(S) = J < d(f),S > dH n+1 L. spt S, 

\S\ = 1 and \d(f)\ < 1 on spt S. Similarly for T — dS. We immediately get that 

|0| = 1, 7i n L spt (T — dS) almost everywhere and \d<f>\ = 1, H n+1 \-sptS almost 

everywhere. □ 

See Figure 4 for an example maximizing form. Note that when S is not unique, this 
proposition implies that \d<f>\ = 1 on the union of supports of all possible minimizing 
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Figure 4: Example of a maximizing form for a square T as input. On the right hand 
side of the figure, details of the form are shown. We visualize the form as a vector 
field. On T — dS, \(p\ = 1 everywhere and on S (red), \d<f>\ = 1 everywhere, while off 
of these sets, both < 1 and \d<p\ < 1. 
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5"s. If we do not have a maximizing form, we have the following modified proposition 
together with the fact that there will be sequence of forms (pi such that T((pi) — >j^oo 
F(T). This easily yields 

Proposition 3.2. Suppose that T(^) -> F(T) = M(S) + M(T - OS). Then there 
is a subsequence of (pi, cp ik such that \(pi k \ — > 1 H n i. spt (T — dS) almost everywhere 
and \d(p\ — > 1 7Y n+1 i_ spt S almost everywhere. 

This modified proposition is necessary since there are easily constructed examples 
having no smooth maximizing form. In fact, the example shown in Figure 4 is not 
actually smooth. The optimizing dcp we show is actually Lipschitz, so it can be 
arbitrarily well approximated by smooth forms even though it is not itself smooth. 
The non-smoothness originates at the points of T — dS where the circular arcs join 
the sides of the square tangentially. At these points, the boundary is merely C 1,1 . 

Remark 3.3. Note that Propositions 3.1 and 3.2 are true if (p and d(p are merely 
measurable wrt the measures 7i n l_ spt (T — OS) and 7i n+l l_ spt S. 

We now state an outline of a reasonable set of steps that would be very useful for 
computations using the dual formulation of the flat norm. 

Conjecture 3.4. Suppose that T is integer density rectifiable with rectifiable boundary 
dT and that both T and dT have finite mass. Define $ to be the collection of all 
Lipschitz forms (p maximizing T((p) and satisfying \(p\ < 1 and \d<p\ < 1. Define X 
to be the closure of the set on which \d(p\ = 1 for every Finally, define § to 

be the collection of all optimizing currents S such that F(T) = M(S) + M(T — dS). 
Then 

(a) There is an S G § such that both S and dS are integer density rectifiable. 

(b) $ ^ 

(c) X = (J spt (S) 

SeS 

(d) \(p\ = 1 on spt(T-dS). 
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Part (a) says that although we are not constraining the minimizing currents to be 
rectifiable, there is at least one in the set of minimizers that is. Part (b) says that 
there are always maximizing forms if we permit them to be merely Lipschitz instead 
of smooth. Part (c) enables us to see where the possible locations for an S might be 
and finally, part (d) gives us T — dS. 

Remark 3.5. We do not expect the refinement and proof of the conjecture to be sim- 
ple. The few previous works in this direction, for example Federer's 1974 paper [5] 
"Real Flat Chains, Cochains and Variational Problems" , are rather technical in na- 
ture. Notice also that this conjecture is only necessary for a rigorous foundation to 
the use of the dual formulation in the computation of the flat norm decomposition. 
Direct optimization over rectifiable currents needs nothing from this conjecture for its 
justification. 

A very simple example where Lipschitz forms are necessary and sufficient for 
optimality is the case in which the current is three equally spaced circles on a torus. 
See Figure 5. Note that the metric on the torus is chosen such that the circles are 
parallel. Figure 6 shows the / of a Lipschitz maximizing form f(x)dy. In the case 
shown of equal spacing between circles, we can't maximize with a smooth form and 
a maximizing sequence must approach the form plotted in Figure 6. In this case, 
the region between the x = and x = a circles or the region between the x = a 
and x = 2a circles is the optimal S. This non-uniqueness is taken into account in 
the above conjecture. Finally, the above propositions, example and conjecture have 
obvious analogs for Fa, the flat norm with scale. 

3.4 L lr TV for co-dimension > 1 

Computing the flat norm decomposition for arbitrary currents permits us to extend 
the L lr TV denoising to sets which are not boundaries or have co-dimension greater 
than 1. One approach is to use the dual formulation of the flat norm. This depends 
on extracting the optimal decomposition from the optimizing form, as discussed in 
the previous subsection. Another approach is to directly optimize over currents. 
We are currently developing both approaches. Figure 7 schematically illustrates the 
decomposition of a 1-current in 3D that results when the flat norm is computed. 
This example actually illustrates both of the useful generalizations possessed by the 
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X 




Figure 5: The 3 circles example. The metric is chosen so that the circles are parallel 
to each other. Optimal S can be either the region between the x = and x = a 
circles or the region between the x = a and x = 2a circles. 




Figure 6: A maximizing form for the 3 circles example. We plot / for the form given 
by f(x)dy. We are forced to use a = 1 and the Lipschitz form plotted here. 
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flat norm decomposition: regularization or denoising of higher co-dimension and non- 
boundary subsets. Notice that the use of the flat norm with scale permits us to choose 
what scale is small and therefore greatly diminished, and what scales are large and 
therefore preserved. In the case of sets which are co-dimension 1 boundaries, we know 
that in a very precise sense, the regularized surface given by spt(T — dS) is the best 
A-curvature approximation to T. See [1 1, 1] for details. 



Figure 7: The green curve is the denoised version of the blue, where we have translated 
the green to make visualization easier. 



3.5 Shape Statistics 

As noted above, the flat norm was previously suggested for shape comparisons in [7, 
13] and then used in [6] for the purpose of computing shape statistics. Our observation 
permits us to use L 1 TV algorithms to compute the flat norm distance for many 
shapes in shape space and the flat norm decomposition that gives this distance. The 
decomposition that we get as a result shows us where the difference is big with respect 
to A and where it is small. See Figure 8. 
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Figure 8: The flat norm via the L lr TV functional provides us with both a distance and 
an informative optimal decomposition into S and T1 — T2 — dS. To use L l TV on sets, 
we need that T\ = dQi and T2 = dQ 2 and either Vt\ C ^2 or ^ C Oi. Using a direct 
method for computing the Fa we don't need these inclusions. Alternatively, we can 
identify the difference between two shapes dQi and dVt 2 as the boundary of the set 
fix A f^2- Then we can use L 1 TV for all codimension 1, boundary shape differences, 
without requiring inclusions. 

4 Summary 

The innovations introduced in this paper simultaneously expand the methods avail- 
able for computing Z^TV minimizers, generalizes L X TV to non-boundary and higher 
codimensional subsets, opens up a new method for multiscale shape decompositions, 
and supports the previous suggestion of Glaunes et al. that the flat norm is useful as 
a distance in shape space. 

Difficulties include the fact that non-boundary O-currents i.e. sets of signed points 
in W 1 not arising as the endpoints of a family of curves, seem rather clumsy to handle. 
For some applications the global curvature bound enforced by the method might be 
too limiting. One can imagine a situation in which the curvature of the approximation 
should be allowed to vary from place to place on the input current or set. One might in 
fact desire something that returns a denoised set whose use as a local mean generates 
a local variance inversely proportional to the locally allowed curvature. 

It is not as easy to handle noise in T in the form of gaps or missing pieces of 
T. Of course, missing pieces means that T is not a boundary. One approach is to 
first denoise dT and then add the resulting S to T to fill in many or all of the gaps. 
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Then T + S can be denoised to remove oscillations by computing the optimal Fa 
decomposition. But then the denoising process requires two steps, each involving a 
choice of A and this seems a bit clumsy. On the other hand this might make sense 
scientifically since the process by which holes and oscillations are generated may be 
quite different with different associated length scales. 

In conclusion the program suggested by the relatively simple observation of the 
relation between L 1 TV and the flat norm promises many new benefits. Many of 
these benefits are immediately accessible while others depend on the some further 
developments outlined above. These, as well as the general expansion of the above 
announcement is the subject of several papers that are in preparation or in planning 
with collaborators. 

A Appendix: Micro-tutorial on Currents and the 
Flat Norm 

If you know something about currents and have a clear picture of the flat norm, 
this section can be skipped. The reference for this section is Frank Morgan's nice 
introduction [9]. The definitive, though formidable, treatise for a fair bit of geometric 
measure theory is still Federer's 1969 tome [4]. References between Morgan and 
Federer include [8, 12]. 

Rectifiable sets Let 7i n denote n-dimensional Hausdorff measure. An n-rectifiable 
subset N of M n+fc is the union of 1) an 7i n negligible set and 2) a countable 
collection of subsets of C 1 n-submanifolds of M n+fc . We are often interested in 
the case where 7i n {N) < oo. Intuitively, an n-rectifiable set looks a great deal 
like an n-manifold at most of its points. 

Currents ra-Currents in R n+k , denoted P n (M n+fc ), are the duals to V n (W n+k ), the 
smooth, compactly supported n- forms on IR n+fc . We will usually suppress the 
]R n+fc and refer simply to T> n and T> n . We restrict ourselves to integer multi- 
plicity rectifiable currents T, which have the following representation: T(<f)) = 
J N 0(ar) ((f>(x),£(x)) dH n ,\/(j) E V n where N is an n-rectifiable set in R n+k , Q(x) 
is an integer multiplicity density function, always ±1 in this paper, <fi is the form 
T is operating on, and £(x) is the unit, simple n- vector defining the orientation 
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on N. Recall that a simple n-vector is the wedge product of n vectors. In our 
case, can be thought of as an oriented representation of the tangent plane 
to iV at x. Changing the sign of the density function has the effect of reversing 
the orientation on N which can also be achieved by replacing £ with — £. 

Currents are naturally equipped with a boundary operator, dT(4>) = T(d4>). 
dT is therefore an (n — l)-current which operates on [n — l)-forms. Note the 
intentional consistency of this definition with Stokes' theorem. 

Intuitively, a current is an oriented manifold or union of oriented manifolds, each 
with an oriented boundary whose orientation is inherited from the orientation 
of the manifold. See Figure 9. There are of course wild beasts in the menagerie 
of currents, but unions of C 1 manifolds with boundary covers a great deal of 
ground, especially when applications are the goal. 




Manifold M, with boundary dM 




The Current Tm showing orientation 
of Tm and OTm 



Figure 9: Example 2-current Tm- Notice that 8Tm = Tqm- The orientation on the 
boundary dM is simply that inherited from the orientation on M. 
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Mass and the Flat Norm The mass of a current is defined as 



M{T) = sup T(4>). (12) 

\(f>\<i,tf>ev n 

Informally, the mass is simply the n-dimensional volume of the rectifiable set 
carrying the n-current. By Theorem 4.1.12 in Federer [4], the flat norm can be 
defined in two equivalent ways: 

F(T) = min (M(S) + M(T - OS)) (given above) (13) 

and 

F(T) = sup (T(4>)) (mentioned above) (14) 

\4>\<i,\dtj>\<i,<t>e'D n 

The corresponding dual definition of the flat norm with scale is given by 

F A (T)= sup (T(<f>)) (15) 
H<i,|#|<A,0e£> n 

Examples of the flat norm decomposition The flat norm involves the optimal 
decomposition of the n-current T into an n-current (T — dS) and an (n + 1)- 
current S. We use the term decomposition in reference to the fact that T — dS 
and S are the components explicitly measured by the flat norm, even though 
T = (T — dS) + dS rather than T = (T — dS) + S. See Figure 10. 

Examples of maximizing forms The dual formulation of the flat norm involves 
finding the supremum over appropriately constrained forms. Figure 11 shows 
a maximizing form for the 2-dimensional disk of radius r. We will discuss the 
computation of optimizing forms and the extraction of S from those optimizing 
forms in section 3.3 

Relation of the flat norm to L 1 The flat norm is very much like the L 1 norm for 
small differences between currents. The value of the flat norm of a difference 
between two currents ends up being roughly the L 1 difference between the close 
parts plus the sum of the n-volumes of what is left. See Figure 12. 
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Figure 10: Example flat norm decomposition. T is the 1-current we are computing 
the flat norm of, and S gives the optimal decomposition into S and T — dS. Finally, 
the flat norm is simply M(S) + M(T — dS) = length of the right-most curve and area 
of the red region. 




Maximizing form (p for disk of radius r 
= fdy — ^dx 



Figure 11: A maximizing form for the disk in 2-D. This form satisfies the constraints 
as long as - < A. The A is, of course, the scale in the flat norm with scale introduced 
above. 
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Figure 12: The flat norm of the difference between these two currents is the sum of 
the area of red region (L 1 like) and the length of the loop that is left. 
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